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An explicit necessary condition for the internal trapping of null geodesies, along with the occur- 
rence of resonance scattering of axial gravitational waves, is proposed for static spherically symmetric 
perfect fluid solutions of Einstein's equations. Some exact inhomogeneous solutions which exhibit 
this trapping are given with special attention to boundary conditions and the physical acceptability 
of the space times. In terms of the tenuity (q = _R/M at the boundary) all the examples given here 
lie in the narrow range 2.8 < a < 2.9. The tenuity can be raised to more interesting values by the 
addition of an envelope without altering the trapping. 
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I. INTRODUCTION 



It is now well known that sufficiently compact stars (polytropic or uniform density) can support the internal trapping 
of null geodesies and the "w-modes" found by Chandrasekhar and Ferrari (1991) ||l|], These modes exist both for 
^ axial and polar perturbations, though the axial ones have been studied more thoroughly. The w-modes in general have 
. no Newtonian counterparts since they are predominantly modes of the spacetime. In the polar case they couple 
■ weakly to the fluid while in the axial case there is no coupling at all. Recent numerical studies of these w-modes have 
^ involved the effect of the equation of state Q and their excitation Q , |6| . Whereas the role that these w-modes may 
(^r^ ' play in real astrophysical processes remains open to much further investigation, it is fair to say that little is actually 
, known about the behavior of the governing potential of the wave equation in exact solutions of Einstein's equations. 
' Such knowledge is important since it is both a route to the physical understanding of relativistic phenomena and a 
check on numerical procedures. The purpose of this paper is to explore necessary conditions for the internal trapping 
of null geodesies and the existence of w-modes (when the centrifugal part of the potential dominates) in physically 
acceptable exact isolated static spherically symmetric perfect fluid solutions of Einstein's equations. We are able to 
exhibit physically acceptable exact solutions which have trapping and which could support w-modes. 



CT; 

. II. REVIEW OF PERFECT FLUIDS 

. ^ ' -^ny metric is an "exact" solution to Einstein's equations. However, the consequent energy-momentum tensor is 

J almost never of any interest. What is of interest are solutions which might have some contact with reality. Recently 0| 

5_] ' a collection of exact isolated static spherically symmetric perfect fluid solutions have been subjected to the following 

I elementary criteria for physical acceptability: 

1. Isotropy of the pressure {p). 

2. Regularity of the origin by way of the scalars polynomial in the Ricmann tensor |^, [p^ . 

3. Positive deflniteness of both p and energy density (p) at the origin. 

4. Isolation by way of the requirement that the pressure reduce to zero at sonic finite boundary radius > 0. 

5. Monotonicity of both p and p to the boundary. 

6. Subluminal adiabatic sound speed (v^ = ^ < 1) 10 • 

Perhaps not surprisingly, only about 10 % of the solutions pass these elementary tests. In what follows we take the 
view that solutions worthy of further consideration must pass all the applicable tests in at least some region |l^ . We 
also take the view that an analytic solution of Einstein's equations can be expected to approximate only a region of a 
realistic configuration. That is, an analytic solution could have an interior causal limit (v^ = 1), a circumstance which 
precludes standard stability arguments [p^ , and yet provide an adequate approximation for a region of a realistic 
configuration. 
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We begin by setting the notation. The hne element in conventional form is (e.g., ||T^) 

= + ^'(de^ + ^HOrdd^') - e'^^^^df (1) 

r 

with the coordinates comoving in the sense that the fluid streamlines are given by u" — e^*^''^i5°. In terms of the 
fimctions <f>(r) and m(r) the regularity conditions reduce to 

$'(0) = m(0) = m'(0) = 0, (2) 

with ' = d/dr and $(0) a constant fixed by the scale of t. Next, in terms of the perfect fluid decomposition 
{T^ = {p{r) + p{r))u°' Ub + p{r)S^) , solving for $'(r) from the r-component of the conservation equations and Einstein's 
equations (VaT" = and G^: - 87rp(r) = 0) we obtain the Tolman ||l6| -Oppenheimer-Volkoff Q (TOY) equation 

= ^P'(^) ^ m(r) + 4:np{r)r^ 

p{r)+p{r) r{r-2m{r)) ' ^ ' 

where, from the t component of the Einstein equations (G^ — — 87rp(r)), 

4M0-^- (4) 

From the TOV equation (here taken to be the right hand members of (^) we observe that p{r) is maximal at r = 0. 
Moreover, if there is an equation of state {p{p)) then either p is maximal wrt r at r = or p is maximal wrt p at r = 0. 
Despite that fact that the TOV equation has been known for over sixty years, only recently has its mathematical 
structure been fully appreciated. For example, we now know that for p{r) > there exists a unique global solution 
for every < p{0) < oo. It is not difficult to find "solutions" of the TOV equation. For example, m(r) can be chosen 
in such a way that (H) yields a solution (with p(r) following from (R)). The simplest choice is clearly m oc but this 
leads us back to the Schwarzschild interior solution. The metric (m) contains two functions, m(r) and ^{r), related 
by (^). The first represents the gravitational energy (effective gravitational mass) (e.g., The second is, in the 

weak field limit r ^ 2m(r), the Newtonian potential. This interpretation offers no insight into the meaning of <I>(r) 
within Einstein's theory, and is a good point to begin our discussion. 

III. NULL GEODESIC LIMIT 

We start with the "centrifugal" part of the potential for non-radial odd parity perturbations. This governs the 
evolution of null geodesies. Radial null geodesies of the metric (0) satisfy 

t = ±[ ^ +D, (5) 

with 9, (f), and D constant. Non-radial null geodesies satisfy 9 — Tr/2 (by choice), 

rV" = 1, (6a) 

0^ 



and 



with 



^V2 = (l (6c) 

r 



B(r)=re-*W, (6d) 
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where • = d/d\ for afRne A, and 6 is a constant > 0, the "impact parameter". The "potential" impact parameter B{r) 
provides, by way of (|6d|) , an invariant physical interpretation of $(r). Null geodesies are restricted by the condition 
&<B(r)g. 

From conditions (||) and the definition (|6d|) it follows that 

B{r) ^ Vor (7) 

as r ^ where V'o is a physically irrelevant scale factor. (The ratio B{r)/b is invariant to scale changes in t.) It 
follows that the necessary and sufficient condition for the internal trapping of null geodesies (that is the existence of 
ro such that r* = and r" < at tq) is given by 

$'(r) > - (8a) 
r 



or, from (||), 



p'^r) < -Bn±m (8b) 



which, with an equation of state (p(/o)) can be given as 

Ar) < _^i!l±ZM . (8c) 
vfr 

From (|8a|) and the TOV equation it follows that 

r < 3m(r) + Anp{r)r^ , (9) 
a relation which makes the trapping of null geodesies a manifestly relativistic phenomenon |23] , p3] . 



IV. FULL POTENTIAL 



The odd parity (axial) w-modes are non-radial perturbations of the spacetime which do not couple to the fluid at 
all. In terms of the frequency zu and mode number I > 2 the governing equation is given by pj] 



where r* is the "tortoise" coordinate 



(^+^ = (10) 



dr^ = — =dr. (11) 

_ 2m(r) 



The potential is conveniently expressed in terms of r and is given by 

- lJ-^2W + 1) + W(p(r) - p(r)) - ^). (12) 
B(rY r 

A necessary condition for the occurrence of resonance scattering of axial gravitational waves by an isolated distribution 
of fluid is a local minimum in V{r) within the boundary of the fluid. (If the centrifugal part of the potential 
( ^(r)^ (?(/ + l)) dominates, which is frequently but not always the case (see below), then (^) provides such a condition.) 
It is the purpose of this paper to explore the occurences of this minimum in physically acceptable exact solutions. It 
is the shape of the function V{r) which is of interest, and since the exterior vacuum has a well known local maximum 
at r ~ 3.28M (for 1 = 2), the boundary conditions associated with the fluid - vacuum interface need careful attention. 
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V. BOUNDARY CONDITIONS 



The Darmois-Israel junction conditions demand the continuity of the first and second fundamental forms at a 
boundary surface. These conditions are well known (e.g., [^) but are usefully reviewed here. We take the "interior" 
metric to be of the form ([^). The "exterior" is the familiar Schwarzschild vacuum (in coordinates (r ^ r, 0, T ^ t)): 

ds^ = + r\de^ + siniOfdcf") - (1 - —)dT\ (13) 

1 r 

At the fluid interface (E), without loss in generality, we take 9 and (p continuous (with intrinsic coordinates 9,(f>,T, 
where t is the proper time). This gives 

rs=r-s. (14) 

The continuity of the first fundamental form is completed by requiring that the particle trajectories at the boundary 
be timelike. The continuity of the angular components of the second fundamental form (extrinsic curvature) give 

M = m(rs), (15) 

and the continuity of the remaining (r - r) component gives 

*E = f ^\.,r, (16) 
rs (rs - 2M) 



which, with the TOV equation, gives 



p(r-s) = 0. (17) 



To summarize, a static spherically symmetric fluid is matched to a vacuum exterior subject to (and only to) (p^), 
(|l|) and (0). Further restrictions are frequently imposed. In particular, if the coordinates are assumed admissible 
(the metric and first derivatives assumed continuous across S) then 

g2*(r.) = 1 _ 2 — , (18) 

and 

ms = = p(rs). (19) 

Whereas ( p^ can be achieved by a simple change in scale (of t or T), in general, (|l9|) does not hold Condition 
( |l8| ) is the necessary and sufficient condition for B to be continuous and continuously differentiable at E. Similarly, 



a simple change in scale makes V continuous but not continuously differentiable at S. The wave equation (10) is of 



course invariant to these changes in scale. In summary, B can be taken to be continuous and continuously differentiable 



at S, and V can be taken to be continuous |27 



VI. EXAMPLES 

Since the uniform density static sphere satisfies (^, one might guess that all static solutions do. This is not the 
case. For example, the Buchdahl solution does not allow a region which satisfies (^. In contrast, the Tolman 
VII solution does (29j . (These are useful exact solutions for the study of the equation of state of neutron stars |Q ) . 
In what follows we demonstrate a number of physically acceptable solutions which do satisfy (^ . We organize the 
examples by way of their motivating ansatz. 
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A. Prescribed form of m(r) 

The Finch-Skea solution is an exact solution which gives reasonable values for the central densities of neutron 
stars. The solution derives from the ansatz 

where C is a constant. The line element can be given in the form 

ds^ = v^dr^ + r^dn^ - A^((C2 - Civ) cos{v) + (Ci + Caw) sm{v)fdt'^, (21) 

where v = VY~\~lu^, uP' = Cr^ and A, Ci and C2 arc constants. Clearly AC2 can be set by the scale of t leaving (say) 
C and /3 = ^ as parameters. The latter is conveniently given by 

p = ~ ^ (22) 



where = -^Z 1 + , or equivalently, in terms of the tenuity a = 

a=-|^. (23) 

The physical restrictions 3 and 6 give, respectively, the following lower and upper bounds to (3 

0.218 <P < 5.605, (24) 

but the limits which follow from B and V are more transparently expressed in terms of a. Up to an irrelevant scale 
factor, the potential impact parameter follows immediately as, 

w 

^^^^ " {l-Pv)cos{v) + {l3 + v)sin{v)' ^^^^ 

We find that B has a local minimum for a < 3 and a local maximum with subluminal sound speed between the local 
maximum and minimum for a >~ 2.768. Some typical plots of B are shown in Fig. |l|. The full potential (up to an 
irrelevant scale factor) is given by 

V{co) = ^ '- ^ (26) 

where 

F( ] = ^ + ^' , 1 {Pv + 1) + (/3 ^ v)tan{v) 
^ ' {(3v ~l)-{(3 + v)tan{v) ' ^ ' 

Some typical plots of V are shown in Fig. ^. We find that there is a local minimum in V (with I — 2) for a <~ 2.933 
and the local minimum lies in a region with subluminal sound speed for a >^ 2.755 ||3^. The Finch-Skea solution 
therefore offers an example of a causal exact solution with trapping |13|. 



B. Prescribed form of ${r) 

A class of models, some of which satisfy conditions 1 through 6, starts with the ansatz 

e2*M =D(l + £;r2)", (28) 

where D and E are constants and n is an integer > 1. The case n = 1 is known as Tolman IV solution |2^]. It follows 
immediately from ( pa| ) that this solution exhibits no trapping {B is monotone increasing and V monotone decreasing). 
For n = 2 condition 3 fails. The case n = 3 satisfies conditions 1 through 6. It has been examined by Heintzmann 
[0, who gives the solution 
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2 l+ar^ 

Again, in terms of the tenuity {a = we find that there is a local minimum in V (with I = 2) for a <^ 2.902 
and the local minimum lies in a region with subluminal sound speed for a >~ 2.788. A typical example is shown in 
Fig. ^ where V has been matched onto the vacuum exterior at the boundary E, and the minimum in V and sound 
speed limit have been indicated. The cases n = 4 and n = 5 satisfy conditions 1 through 6 and have been solved by 
Durgapal js^] following the formulation of Korkina . For n — 4 the solution is given by 



ds' = ^ ,_,,Jii^J ^ '^\cr^ , + ~ + Cr^ fdt\ (30) 

( 7 1 2- ) 



In this case we find that there is a local minimum in V (with I = 2) for a <^ 2.892 and the local minimum lies in a 
region with subluminal sound speed for a >~ 2.780. Since B(r) has a local maximum up to a < 3, it is clear that 
the dynamical part of the potential V can dominate. For n = 5 the solution is 

y 112 ~^ ^l+6Cr^' 

We find similar results in this case. There is a local minimum in V (with I =^ 2) for a <^ 2.886 and the local minimum 
lies in a region with subluminal sound speed for a >^ 2.776. The solutions (p9|), (|o|) and ( |3l] ) cannot represent the 
core region where, we note, they are acausal |l3|. 



VII. DISCUSSION 



Condition (|^) is proposed as the necessary condition for the internal trapping of null geodesies, and for the occurrence 
of resonance scattering of axial gravitational waves when the centrifugal term dominates the potential, in static 
spherically symmetric perfect fluids ||3^ ]. This condition is not always satisfied. For example, it is not satisfied in the 
Buchdahl 128] solution. We have demonstrated some physically acceptable exact solutions for which the condition is 
satisfied. One, the Finch-Skea solution, offers an example of a complete causal exact solution with trapping. At the 
very least, these examples can provide a check on numerical procedures which attempt to gauge the role that w-modes 
may play in real astrophysical processes. In every case studied here we have found that for resonance scattering the 
tenuity {a = r^/M) lies in the small range 2.8 < a < 2.9. (In all cases, as a decreases, the causal boundary (v^ = 1), 
if it exists, moves out and approaches the minimum in V for some minimum a, exactly as expected.) Whereas this 
range is well above the Buchdahl limit of 9/4 |3^], it is too low for, say, neutron stars (Unphysical solutions with 
trapping and a > 5 are known [^9| .) It is reasonable to suggest, as has often been done, that an exact solution may 
reflect only part of a more realistic configuration. Boundary conditions within a distribution are easily derived from 
the discussion given in section V: p{r) and m(r) must be continuous to avoid surface layers (shells). In particular, 
m(r) need not be continuously differentiable (which, at least formally, allows first-order phase transitions). All that 
is needed to raise a into a more interesting range (say 3 < a < 10) is the addition of an envelope. The envelope is 
constructed subject to the continuity of p{r) and m(r) at Va-, where r^r is exterior to the minimum in V, and must 
allow p(r^) = at a finite boundary where < rs | |4C| ]. 
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FIG. 1. The potential impact parameter B for the Finch - Skea solution. The curves (bottom to top) have {a,WY;,l3) given 
by (3,^,2.638), (2.856, 1.529,4) and (2.768, 1.614,6.332). The vertical scale is irrelevant. For a < 2.768 the sound speed is 
superluminal beyond the local maximum in B. The curves are continued to E. 



FIG. 2. The potential V (for mode number I = 2) for the Finch - Skea solution. The curves (bottom to top) have (a, we, /?) 
given by (2.933, 1.464,3.105), (2.856, 1.529,4) and (2.806, 1.575,5). The vertical scale is irrelevant. The curves are continued 
to S. 



FIG. 3. The potential V (for mode number I = 2) for the Heintzmann solution. Here we use a = jj. as = 2.850. V has a 
local minimum at a = 2.717 and the sound speed becomes superluminal below a = 2.291. V for the vacuum exterior is also 
shown. 



8 



This figure "fsb.png" is available in "png" format from: 



http://arXiv.org/ps/gr-qc/0007073v3 



This figure "fsv.png" is available in "png" format from: 



http://arXiv.org/ps/gr-qc/0007073v3 



This figure "heint.png" is available in "png" format from: 



http://arXiv.org/ps/gr-qc/0007073v3 



